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Abstract 

In the Hamiltonian language we provide a study of flat-space cos¬ 
mology in bigravity and massive gravity constructed mostly with de 
Rham, Gabadadze, Tolley (dRGT) potential. It is demonstrated that 
the Hamiltonian methods are powerful not only in proving the absence 
of the Boulware-Deser ghost, but also in solving other problems. The 
purpose of this work is to give an introduction both to the Hamilto¬ 
nian formalism and to the cosmology of bigravity. We sketch three 
roads to the Hamiltonian of bigravity with the dRGT potential: the 
metric, the tetrad and the minisuperspace approaches. 


*E-mail: Vladimir. Soloviev@ihep.ru 


1 



Contents 

1 Introduction 

2 Metric variables 

3 Tetrad variables 

4 Minisuperspace 

4.1 Gravitational potential terms . 

4.2 Gravitational kinetic terms. 

4.3 Matter terms. 

4.4 List of GR formnlas given for reference. 

5 Matter couplings 

5.1 f-matter and g-matter . 

5.1.1 Massive gravity on Minkowskian backgronnd . . . . 

5.1.2 First branch of bigravity. 

5.1.3 Second branch of bigravity. 

5.2 g-matter withont f-matter. 

5.3 One matter minimally interacting with two metrics . . . . 

5.4 One matter minimally interacting with the effective metric 

5.4.1 First branch. 

5.4.2 Second branch. 

6 Non-dRGT bigravity 

6.1 Non-dRGT massive gravity . 

7 Conclusion 

8 Appendix 


a 

s 

16 

12 

18 

19 


20 

23 

23 

25 

26 

27 

28 
32 
34 


p 

Ei 


37 


37 


1 Introduction 

One of the actual problems in modern physics is the problem of dark energy. 
There are different approaches to it, and one of them is to revise the standard 
theory of gravitation, i.e. the General Relativity (GR). New predictions are 
expected for physics at the very large distances, and they may be provided 
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by making gravity massive. The popular way to do so is to introduce a new 
tensor field with the same properties as the metric. From the aesthetic point 
of view it is natural to expect that this new field will be dynamical. Such a 
model is called bigravity or bimetric gravity. 

The crucial problem in constructing the bigravity or massive gravity mod¬ 
els is to avoid the Boulware-Deser ghost pQ. This problem was solved by de 
Rham, Gabadadze, and Toley |21I3] (dRGT), who proposed the genuine form 
of potential involving the matrix square root. The proof that this theory 
is ghost-free has been given first in the metric Hamiltonian approach mu- 
But it is hard to work with the dRGT potential in metric formalism, whereas 
in tetrad approach the calculations become explicit. Fortunately, the back¬ 
ground cosmology that will be discussed in this work exploits only rather 
simple form for both the two metric tensors, and in the so-called minisuper¬ 
space approach the relevant matrix is diagonal, so the square root of it can 
be found directly. 

Let us mention that the Boulware-Deser ghost appears below even in 
minisuperspace models for the two cases: first, when a non-dRGT potential 
is chosen (see Section 6), second, when one matter field minimally couples 
to a pair of metrics (see subsection 5.2). 

A lot of work in cosmology of bigravity has been already done by many 
authors 0. As a rule, the main instruments were the Lagrangian equations 
and Bianchi identities. The aim of this work is to demonstrate the main 
features and power of the Hamiltonian formalism in analyzing the variety of 
cosmological problems from a unique viewpoint. This variety arises due to 
many ways by which the matter fields can couple to the gravitational fields. 
We consider the minimal interactions of two species of matter with the two 
metrics, the double interactions of one matter with the two metrics, and the 
effective metric interacting with the matter. We also illustrate the general 
scheme by an example of the non-dRGT potential. For simplicity we consider 
only fiat space dynamics, and each metric has only one matter source. 

Here all the local physics is ignored including the important question of 
stability of the background cosmological solutions. This is discussed in many 
articles, for example, see recent publications PE]. We are impressed by the 
statement [8] that it is possible to cure the growing perturbations by taking 
a tiny Planck mass for the second metric 

^We are sorry for citing only part of them and not trying to give a complete list of 
references. 
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In Sections 2 and 3 we briefly explicate the constrnction of the bigravity 
Hamiltonian (first, in metric variables; second, in tetrad variables), for a 
technically more detailed presentation we address the reader to articles |9l 
For approaches to the Hamiltonian formalism that are not 
discussed here see H El ng. The title of this work should remind about 
a good old book na. The mini-superspace Hamiltonians are constructed in 
Section 4, the machinery given there is applied to the concrete derivations of 
the cosmological equations in Sections 5, 6. 


2 Metric variables 


The Lagrangian density of the bigravity theory may be written as follows 

‘ 2771 ^ 

- U 

K 

where the two copies of the GR Lagrangian are involved 

The matter Lagrangian and the potential for interaction of the two metrics 
are also included in C. For simplicity, as an example of matter, we will take 
the scalar held with a minimal interaction with some metric 

^ . ( 2 . 1 ) 


The role of metric below may be played by or by or by their combina¬ 
tion called the effective metric (see subsection 5.3). The constant k standing 
in front of the potential may be chosen coinciding with standard Newtonian 
one Kg = IGtiG, but in some works it is different, for example. 


( 2 . 2 ) 


1 _ 1 1 

K Kg Kf^ 

SO we start with the general notation. We denote 3-1-1 components of a 
general metric and its inverse in the ADM coordinate basis as follows 




( -AA-2 \ 

V ) ’ 


(2.3) 

(2.4) 
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where M is the lapse, TV* is the shift vector, -0^ is the induced metric on 
spatial hypersurfaces of fixed time, and is fhe inverse matrix for it. The 
interaction potential of the two metrics suggested by de Rham, Gabadadze 
and Tolley [21 13] is constructed as a linear combination of the symmetric 

polynomials e* of matrix = (y\/g~^f^ 

eo = 1, 

Cl = Ai + A2 + A3 + A4, 

62 = A1A2 + A2A3 + A3A4 + A4A1 + A1A3 + A2A4, 

63 = A1A2A3 + A2A3A4 + A1A3A4 + A1A2A4, 

64 = A 1 A 2 A 3 A 4 , (2.5) 


where A* are eigenvalues of X. Then 

4 

U = /3„e„(X) = + ... + —f = NU , 

n=0 

Our plan of constructing the Hamiltonian formalism is the following. 
First, we import the canonical variables and the expressions of Hamiltonians 
from the two copies of General Relativity 


Hg = j d^x {NH + N^Hi) , 
Hf = j {NU + N^Ui) , 

where 

= - 24 , 

Hi = -2m.. 

* A] 


( 2 . 6 ) 

( 2 . 7 ) 

( 2 . 8 ) 
( 2 . 9 ) 

( 2 . 10 ) 

( 2 . 11 ) 


The canonical gravitational variables have the standard Poisson brackets 


{7.j(x),7r“(j,)} = ^-{6}6‘ + 6‘6f)6{x.y), (2.12) 

{ry,{x),n“(y)} = ^-{SfS‘ + S‘Sl]S{x,y). ( 2 . 13 ) 
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The ADM-decompositions for the two space-time metrics are as follows 


II/..II = 
iini = 

\\9n\ = 

Next we construct the matter Hamiltonian 


+ rimnN^N^ VjkN’^ \ 
VikN^ Vij ) 

_iV-2 ATiAT-s \ 

iVW-2 - N^Nm-^ J ’ 

-iV2 + 7^_„iV-iV- 

likN^ lij ) 

_Ar-2 


iVW-2 7*i - ) ’ 


Hm = / d^x (^^nM + 


(2.14) 

(2.15) 

(2.16) 
(2.17) 


and in case the role of matter is played by one scalar held with a minimal 
coupling to one metric according to Eqs. (12.31) . (12.41) we obtain 

+ \/^W(0). 

In dealing with two or even more metrics instead of the coordinate basis it is 
suitable to use the geometrical basis [HI UHl [201 |2l] formed by a unit normal 
to the state hypersurface n°‘ and three vectors tangential to this hypersurface 
ef. Let be an arbitrary space-time coordinate system, and (t, x*) is the 
ADM coordinate system, i.e. functions = X^{t) give one-parametrical 
foliation of the space-time by spacelike hypersurfaces, whereas x® are internal 
coordinates on a hypersurface. Then ef = dX^/dx^ are the three tangential 
vectors. Of course, we need a space-time metric to construct the unit normal, 
let us suppose this metric to be then 

= -1, f^,un^ei = 0. 


If we decompose over this basis another metric, for example, we obtain 


11^"^ 

I 


f —u ‘^[n^n''] u^u \ 

/ {-u^ +'jmnuV)[nf,n^] \ 


(2.18) 

(2.19) 
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where new variables m, appear. Their role is, at least, three-fold: first, 
they appear in the relations between two pairs of lapse and shift functions: 


u = 


N 

iV’ 


= 


N^-N^ 

N 


( 2 . 20 ) 


second, they are formed by projections of tensor onto the basis (ua, 
constructed with the metric 


u = 


u = 






( 2 . 21 ) 


third, they are the coefficients that connect the two bases (ha, e),,) and (na, e),): 


u 


n,, = un,,, e;, = e:, — = —n^ — —ef, ef = ef. ( 2 . 22 ) 


''M) 


1 

u u 

As the bases (n^, e)^) and (n“, e“) are formed with the help of metric 
this metric has only 6 nontrivial components in decompositions over them 


\\r 


/ —[n^rf] 

V 0 

( -[nf,nj\ 
\ 0 





(2.23) 

(2.24) 


One can see that all the components of matrix = g^^fau are functions 
of variables u, m*, 7 *^ and gij, the same statement is valid for the invariants 
of this matrix, for example, Tr(Y"). Therefore the eigenvalues of X = \/Y 
are also dependent only on these variables. In general case it is impossible 
to obtain the explicit expression for the dRGT potential in the form 


U = NU{u,u\gij,jij). 

For this reason we start with a general function U and work with its formal 
derivatives 


V 

V. 


dU 

dU 

(9m*’ 


(2.25) 

(2.26) 
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together with the function 


~ dU ,dU 

The Hamiltonian now appears in the following form 


H = / dP'x 


2m^ 


N in + un + + nM + [Ui + Hi + UMi) 


Then the derivatives of the total Hamiltonian over non-dynamical variables 
N, N^, u, n* provide us with the primary relations: 


7^ = 
7^. = 
5 = 

5. = 


OH 2??7,^ 

tttt = H -\- uH + HHi + Hm H- U = 0, 

ON K 


OH 

1 OH 

N~du 

1 OH 


= Hi + Hi + Hmii 

OHm 2m^ 


= H 


du 


+ 




K 

2m? 


H = 0, 


N du^ du'' u 

Then we can rewrite the Hamiltonian as follows 


= 0 . 


H = d^x [N (7^' + + u%) + N^TZi] , 


where 


2t?7^ 

n' = H + - — w, 


K 


For the potential u of the general form equations 

iS = 0, Si = 0, 


(2.27) 

(2.28) 

(2.29) 

(2.30) 

(2.31) 

(2.32) 

(2.33) 


may be solved for variables u, u\ and then these solutions should be sub¬ 
stituted into the expression of 71' in order to make it a constraint, i.e. a 
function of canonical variables only. Ater that equations 


7 ^ = 0 , = 0 , 


(2.34) 


occur 1st class constraints arising as a consequence of the space-time coor¬ 
dinate invariance of the bigravity. On-shell, i.e. on the surface of variables 
satisfying Eqs. fl2.33p they form the celebrated Poisson bracket algebra 

{7^(x), 7^(|/)} ^ {r,'\x)nk{x) + 5,(x, y), (2.35) 
















{Tii{x),'Rk{y)} ^ Tii{y)5^k{,x, y) + '}Zk{x)6^i{x, y), 


(2.36) 


7^(a;)5^(a;, y), (2.37) 

discussed by Dirac in the proposal to construct the held quantization on 
the curved surfaces [22]. For the potential U of the general form we obtain 
equations necessary to fulhl the above Poisson algebra 


dr]jk 

n i 9U , dU 
d'yik du 


o du ^ dii 
“ 1 “ 


djjk 


+fr 

u-' 7 — - bjU 
du^ 

= 0, 

(2.38) 

i dU 

-“A 

= 0. 

(2.39) 


In order equation iS = 0 to become a geniune constraint excluding the 
Boulware-Deser ghost degree of freedom it is necessary to put one more con¬ 
dition on the potential. This condition is the requirement that it should be 
impossible to solve Eqs. (l2.33p for auxiliary variables = (m,m*), i.e. these 
equations are to be functionally dependent, 


D{S,S,) 

D{u, M*) 


det 


du^-dvf’ 



u + n 


M 


= Hess(4x4) 



U + n 


M 


0 . 


(2.40) 

Next we need to prove the existence of a secondary constraint D, and 
the fact that it has a non-zero Poisson bracket with S. For this purpose it is 
useful to apply a method suggested by Fairlie and Leznov [23] on constructing 
implicit solutions of the homogeneous Monge-Ampere equation. This trick 
has been hrst developed and applied to massive gravity by Comelli, Nesti, 
and Pilo [H]. Combining Eqs. fl2.38p . fl2.39p with the technique of Ref. [23] it 
is possible to prove that 


{5(a;),5(2/)} 0, 


and therefore 


5 = {5,H} 


d^xN{S,n'} 


d^xNQ = 0 -f-)- D = 0. 


By means of the Jacobi identity it is also possible to show that in general 


{D,5} 7 ^ 0 , 
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and so these two constraints are second class in the Dirac terminology. The 
anxiliary variable u is hxed by the condition of preservation the secondary 
constraint in time 


{D,H} 


d^xN ({D, 5}) = 0. 


The other anxiliary variables n* are to be determined from eqnations Si = 0. 
For this pnrpose it is necessary to reqnire 


det 


dS,, 


du^ 


= det 




du'^du^ \ K 


2m^ 


u + n 


M 


= Hess 


(3x3) 


2mr 


K 


U + Hm 7^ 0, 


( 2 . 41 ) 

and therefore the rank of the big 4x4 Hessian fl2.40p shonld be eqnal to three. 
Finally, we come to the Hamiltonian formalism of bigravity with 12 pairs of 
canonical variables fonr 1st class constraints 7l,7li, and 

two 2nd class constraints S, D. There are 7 gravitational degrees of freedom 
for snch a gravitational potential. 

Here we followed articles P[IQ]. For the massive gravity case see also m- 


3 Tetrad variables 

It is possible also to constrnct the Hamiltonian approach to bigravity with 
the tetrad variables connected with the metrics by formnlas 

S"- = ESBVi"'®, = Fi^FfhAB, (3.1) 

where Hab is the Minkowski metric 

Hab = diag(-l, 1,1,1), 

and (i7^, E^), F^) are the mntnally inverse tetrad matrices. In general 

the tetrad Hamiltonian formalism is more involved than the metric one. As 
10 components of the metric tensor are replaced by 16 tetrad components the 
nnmber of constraint eqnations increases. Therefore the constraint algebra 
becomes more extended. In the same time the tetrad representation allows 
to obtain an explicit expression for the dRGT potential. Indeed, the matrix 

= E^^^F^a, (3.2) 
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occurs a square root of the matrix = g^°‘fav, if the symmetry equations 

E>XF;^ - = 0, (3.3) 

are valid. Then instead of Eqs. (12.51) it is suitable to apply equivalent formulas 
for the symmetric polynomials 


Co 

ei 

62 

63 

64 


1 , 

TrX, 

i ((T1X)2 - TlX"), 

i ((TiXl^" - STrXTrX^ + 2TrX®) , 

6 

^ ((TrX)^ - 6(TrX)2TrX2 + 3(TrX2)2 + STrXTrX^ 
detX = 


det IlF^all ^ 
det||E^a|| v^' 


6TrX^) = 

(3.4) 


It is convinient to start with the privileged form of a tetrad, and a bit later 
transform it to the general form. Starting hrst with the metric we pro¬ 
pose to use as such a simplihed tetrad the following set of four space-time 
vectors: the first vector Eq is the unit normal to the hypersurface of state 
ft" constructed on the base of metric g^i,: 


g^uU^^e^ = 0 , 


g/iun^rf = - 1 , 


the other three space-time vectors are the three tangential to the hypersurface 
vectors provided by the triads of the induced metric y^j-For this triads e“ we 
have equations similar to Eqs. fl3.ip : 

'Jij ^i^j^abi ^ab diag(l, 1, 1), ^j' 


Then the tangential space-time vectors are formed as follows 


E" = eV", 


e" = e" 


ax" 

dx^ 


The inverse tetrads are 


El = 


El = 
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Following this way for the metric we also get a priviledged tetrad 

J-” = -n„ (3.6) 

4 = (3'6) 

The general tetrad arises as a result of a boost transformation 

= ( pa pa ^ ^ = Vl+PaP“, K = > (3-7) 

of the priviledged one = A^^{p)F^. Here 

p“ = Pa, = hcD- 


Then we can derive the matrix dehned in fl3.2p . and obtain the following 

_ / A[n^ny\ 


= 




'1^3 i 


(3.8) 


where 


A 

Ci 


e 

u 

Pa^ 

U 

Fe 


u 


-PaG 


M*pafic 


U 


+ F“Pafee 


ib 


(3.9) 

(3.10) 

(3.11) 

(3.12) 


After straightforward calculation of all the symmetric polynomials fl3.4p we 
get an explicit expression for the dRGT potential which occurs linear in 
variables u, F 

U = uV + u% + W. (3.13) 

Functions V, Vi, and W are given in Appendix. They depend on canonical 
variables eia, ha and on auxiliary variable pa- Especially, V) occurs linear in 
Pa 

Vi = -PaCabfib- (3.14) 

As the potential term does not contain any velocities, it enters the Hamilto¬ 
nian of bigravity without any changes. 
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Next we are to find the corresponding Hamiltonian form for Cg and Cj. 
This is familiar from GR [21]. The new canonical variables are triads ej^, 
together with their conjugate momenta 7r®“, n*“, and their Poisson brackets 
are the following 

{fia(a;),n^'’(2/)} = {eia{x),n^^{y)} = 5l5l5{x,y), (3.15) 

{^ia{x),eji,{y)} = {7r“(a;),n^'’(2/)} = {fi„(a;),f,b(|/)} = 0, (3.16) 

{eia{x),ejb{y)} = W{x),n^\y)} = {Yi"^{x),IV\y)} = (3.17) 

Momenta of the metric formalism are expressed through triad momenta as 
follows 

nP = 1 + F“n“), (3. is) 

TT*^' = ^ (e*V“ + eV) . (3.19) 

Then the Poisson brackets of the metric momenta are nonzero off-shell 

{n*^(a;), Ti^^{y)] = ^ , (3.20) 

{n^^{x),n^^{y)} = ^ •(3-21) 

Here the new constraints specihc for the tetrad approach appear 


M” = = 1 (f‘°ni - F“ny = o, (3.22) 

= iUte’V* = 1 (e“< - <JX) = 0, (3.23) 

or in a different form, 

Lab = f^a^l~f,bK = 0, (3.24) 

J^ab = eiaY - eibK = 0. (3.25) 


The combined Hamiltonian, of course, contain these constraints 

Bg^f = Bg + Bf = 

= j <Yx {NH + N^Hi + r^Lab) + 

+ [ Yx {NU + WUi + X^^Lab) , (3.26) 
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together with the constraints Ti, Hi, Ti, Hi known from the metric Hamilto¬ 
nian fl2.8p - fl2.1ip . As in bigravity the two metrics are mixed in the poten¬ 
tial, only diagonal rotations of triads Cja, leave the Hamiltonian invariant. 
Therefore symmetric combinations 


— Lab + Lab — 0, (3.27) 

become 1st class constraints, whereas constraints 

^ab ~ “ Lab = 0, (3.28) 


are 2nd class. 

Finally, we are to take into account the symmetry conditions 
the Hamiltonian variables they take the following form 

Ga = Pa + - U^Vabij = 0, 

Gab = ^icPcla^b] — O' 

The Hamiltonian of bigravity now appears in the following form 


H = 


HxN 




K 


K 


K 


+ / Sx [N\l-Li + Hi) + -F Xab^ab + L'^Ga + A“'’G'ab] . 


fCT . In 

(3.29) 

(3.30) 


2 

-w 


+ 


(3.31) 


When two types of matter, f-matter and g-matter, separately and minimally 
couple to metrics and their contributions can be incorporated into 
constraints H, Hi, H, Hi. Besides the canonical variables the Hamilto¬ 
nian depend on the Lagrangian multipliers A“, A“^, providing Eqs. 

fl3.27p - fl3.30p and on other auxiliary variables. Varying over u, u^, N, V® we 
obtain the following equations 


5 

n 

TZi 



2m\^ 

H + 

V : 


K 



Hi + 

- 


K 


TZ' + uS + vZSi = 0, where 7Z' 
Hi -|- Hi = 0 . 


2777 ^ 

H + ^W, 

K 


(3.32) 

(3.33) 
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(3.34) 

(3.35) 


















We do not vary over variable Pa that enter fl3.3ip in non-linear way. Eqs. fl3.29p 
can be solved for m® 



(3.36) 


where 



(3.37) 


Eqs. fl3.33p contain variable Pa in a linear way and so can be easily solved 
for it 



(3.38) 


Eqs. fl3.30p depend on the canonical variables e^, he, and on the auxiliary 
variable Pa- After the exclusion of Pa Gab become constraints. They have 
nonzero Poisson brackets with and together form 6 constraints of the 
2nd class. 

In the tetrad Hamiltonian formalism there are 18 pairs of the canonical 
variables (eja,7r®“), (fja,n®“), 7 constraints of the 1st class 77', 77*, and 
8 constraints of the 2nd class S, H, Gat- Of course, we again have 7 
gravitational degrees of freedom. The auxiliary variables pa are found from 
equations iS* = 0 which are linear in pa- Variables m® are found from equations 
Ga = 0 linear in m®. The last auxiliary variable u can be found from 



(3.39) 


this is also a linear equation. 

The above exposition followed articles da [13]. Different treatments can 
be found in Refs. [HI ESI ESj- 

4 Minisuperspace 

Let us look at the cosmological problem, i.e. consider the Hamiltonian evo¬ 
lution of a homogeneous and isotropic Universe in the bigravity theory. For 
simplicity and brevity only the spatially flat case will be treated. Most of 
attention will be paid to the theory with the dRGT potential. 

In general, the bigravity Hamiltonian is a sum of kinetic and potential 
gravitational terms, composed of two metrics and plus one, or two 
matter contributions. 
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We take the Friedmann-Robertson-Walker ansatze for the two metrics 


f, u = = (4.1) 

g, . = v^=Nue, ( 4 . 2 ) 

and introduce a new variable r = u/^. 


4.1 Gravitational potential terms 

For the admitted metric tensors of the form fl4.ip . fl4.2p we obtain 

9~^f = g^'^fau = diag . (4.3) 

Let us suppose 

> 0, M > 0, ^ > 0, o; > 0. (4.4) 

In fact, in calculating the matrix square root of the positive diagonal matrix 
fl4.3p we consider only one (positive) square root for each positive expression [§ 

X = \/g~^f = diag (^+\/u-2, +\/r^(5y j = diag r6ij) , (4.5) 

and eigenvalues of this matrix are determined from the following equation 

det(X - AJ) = 0. (4.6) 

It is easy to hnd that 

Ai = M A 2 = A 3 = A 4 = r, (4-7) 

and then construct the symmetric polynomials: 

eo = 1 (4.8) 

Cl = Ai + A 2 + A 3 + A 4 = - + 3r, (4.9) 

u 

3?^ 

62 = A 1 A 2 + A 1 A 3 + A 1 A 4 + A 2 A 3 + A 2 A 4 + A 3 A 4 =-h 3r^, (4.10) 

u 

3r^ 

63 = A 1 A 2 A 3 + A 2 A 3 A 4 + A 1 A 3 A 4 + A 1 A 2 A 4 =-h (4.11) 

u 

3 

64 = A 1 A 2 A 3 A 4 = —. (4.12) 

u 

^In this special case nothing serious will change if we take another sign for some eigen¬ 
values. 
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The dRGT potential is a linear combination of these polynomials which oc¬ 
curs a linear function of variable u: 

4 4 

Pid = Nu^ = N{uV + W), (4.13) 

where 

R = + + + = (4.14) 

W = /3if + 3/32^^u; + 3Mu;^ + /34u;^ = ^^Bi(r), (4.15) 

Bi(r) = /3i + 3/3i+ir + 3/3i+2r^+ /3i+3r^. (4.16) 


4.2 Gravitational kinetic terms 

If we start with the Hilbert-Einstein gravitational Lagrangian density for the 
space-time metric 

Cg = (4.17) 

Hg 

and exploit the cosmological ansatze 




(4.18) 


we get 


Rg — Tg — — 




6aW^2 


(4.19) 


where 77 is the Hubble constant. 

Dehne the momentum canonically conjugate to a 

dCc 12a. 12a2 niia r i i 

It is suitable to use H instead of Tia, then we will use non-canonical Poisson 
bracket: 


{a,H} = 


'l2a2‘ 


(4.21) 


In bigravity we get the kinetic part of the Hamiltonian in the following form 


^kinetic Tf 7g 


Gu^N 2 G^^Nu 2 


(4.22) 
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This part generates the kinematical Hamiltonian equations 

^ ^kinetic} NujHj^, ^ Hkinetic} (4.23) 

equivalent to the definition of two Hubble constants Hf, Hg. It is easy to 
derive an equation for the evolution of the relative scale factor r 

r = Nr{Hf -uHg). (4.24) 


4.3 Matter terms 

As an example of matter, we take the scalar field having a minimal interaction 
with gravity 

- W(0)^ , (4.25) 

where is the relevant metric of the form given by Eq. (I4.18p (there will be 
different cases below), where 0,^ = dcp/dx^. In the homogeneous cosmology 
0 = 0(t) and 

= 70[tt^] + U^. (4.26) 

Scalar fields may be of special interest in studying inflation in the bigrav¬ 
ity m If this matter field interacts with only one metric there is a 
straightforward correspondence between the velocity and the momentum 


Ts = 


W’ 




d% 

00 


Jf"' 


M 


(4.27) 


We can also introduce density p and pressure p of matter according to the 
formulas 


= Afa^p, 




(4.28) 


Below we will also consider the case when a scalar field simultaneously couples 
to metrics and fgiy. Then 


70 


^^0^ a;^0^ 

~W ^ 2N ’ 




Ntta 


^3 n ’ 


(4.29) 
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and therefore 


Ts = 


Nnl 


2 ( ^ + wsV 


U^ = N + 0 ;=^) W, 


(4.30) 


= N{ue + 


T^a 


2 ( S4 + Cj3 

* U 


3 + a;3) 


+ W(0) 


Here T-La is a nonlinear function of u. 


(4.31) 


4.4 List of GR formulas given for reference 

^9 = —V^{9^''R9U - 2 A) + >Cm( 0 ^, g^u), K-g = IOttG, 

Kg 

^<P = - W( 0 )^ , 

- -g>^'^(R - 2A) = _ 

2 2 dggu 

g,u = {-N^{t),emj), =Ne, 0 = m- 


(4.32) 

(4.33) 

(4.34) 

(4.35) 


^9 = Ne - 


3 1 6 / e 


live 


- 2A 


^ „ 
Hq = , TTf =- Hg. 

^ A^e’ ^ Kg ^ 



i! 

at' 


7r0^ 

2^6 




/vn 


12^2 


(4.36) 

(4.37) 

(4.38) 

(4.39) 
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2 


(4.40) 



2A 

+ P+ — 


K 


9 


0 , 




SnG 




i = mHg, (4.41) 

Nk 

Hg = - ^{p + p) =-4:7iGN{p + p), (4.42) 

p = -3^{p + p). (4.43) 


If A and the matter equation of state p = p{p) are known, the initial condition 
can be taken for p or for Hg. N is arbitrary in accordance to the freedom of 
time reparametrization. 


5 Matter couplings 

5.1 f-matter and g-matter 

The hrst works on cosmology of bigravity [2H1129], [301 El] appeared rather 
soon after the proposal of the bigravity theory with the dRGT potential made 
by Hassan and Rosen |1]. In trying to solve the dark energy problem it is 
enough to admit g-matter only, but if also dark matter is under study then 
f-matter is also required mi. 

If there are two species of matter: f-matter and g-matter, and each one 


has a minimal coupling to or then 



H ^^potential ^^kinetic ^^matteri 

(5.1) 

where 



^^potential 

N{uV + W), 

K 

(5.2) 

^^kinetic 

= —N + —Nu , 

K.f ■’ Kg ^ 

(5.3) 

Hmatier 

= Noj^pf + Nu^^pg. 

(5.4) 

The Hamiltonian can 

also be written as follows 



H = A^7^ = NW + NuS, 

(5.5) 
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where primary constraints are 1Z = TI' + uS and 5, or better, VJ and S: 


n' = = 0, (5.6) 

Kf ■’ K 

S = + + (5.7) 

hig hi 

For the scalar sonrces we have 


Pf - 2cn6 

(5.8) 

Pg = 

(5.9) 


The constraints fl5.6l) . fl5.7p may be written as a conple of Friedmann eqna- 
tions: 




Ht = 


+ T' 


Af = m 


2^/ 5i(r) 


K 




-pg + 




A = rn^—Bn(r) 

K 


(5.10) 

(5.11) 


We see that, in contrast to GR fl4.4Up . the cosmological terms natnrally ap¬ 
pear in the bigravity, and they are not constants, bnt dynamical qnantities. 
Hamiltonian eqnations for the matter helds provide ns with a conple of con¬ 
servation laws 


Pf + ^-ipf+Pf) = 

OJ 

p 

Pg A ^-^{Pg + Pg) = 0. 


(5.12) 

(5.13) 


The secondary constraint arises in order to preserve the primary con¬ 
straint dnring the evolntion 


S = {S,ll} = N{S,n'} = NQ = 0, (5.14) 


where 

f)T?7 ^ 

G = {5, n'} = - {uHf - iHg) -t- = 0. (5.15) 

hi 
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As the secondary constraint is factorized 

^ 2 5 

there are two branches of solutions for it: 

= 0, Hg = tHj, 

and 

= 0, ■(-)■ /3i + 2/?2'r + = 0. 

Below we consider both two cases. 

Variable u is hxed by the following requirement 

n = {n,R} = N ({12,7^'} + m {12, 5}) = 0, 

here we also have two cases: 


(5.16) 

(5.17) 

(5.18) 

(5.19) 


u = 




or u = 


{122,7^'} 


( 111 , 5 }’ ^ { 122 , 5 } 

For the two Hubble constants we get dynamical Hamiltonian equations 


(5.20) 


Hj = {Hf,}i} = {Hf,u} — , 

(5.21) 


(5.22) 


which appear in the following explicit form (with account for Eqs. fl5.6p . fl5.7p i 


Hf = 


H, = 


Nk 


f 


4 

Nuk, 


Pf +Pf + {I— ur) 
Pg+Pg- i^-ur) 


2nn? Di{r) 

K 

2rn? Di{r) 


K u 


where 


Di{r) = A + 2A+H + A+ 2 r^ 


(5.23) 

(5.24) 

(5.25) 
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5.1.1 Massive gravity on Minkowskian background 

If we suppose to be a non-dynamical Minkowskian metric 

N = l, u=l, r = ^, (5.26) 

and, of course, exclude f-matter 


Pf = 0, Pf = 0, 


then we get a massive gravity theory 


H 

5 

Q 


uS 

e 


2m? 


K 


-W, 


6 2 n A 

+ —Bo 


12 m^ 


K 


= 0 , 

^H,(/3ie + 2/32^ + /33) =0. 


(5.27) 


(5.28) 

(5.29) 

(5.30) 


From the last equation it follows that there are no non-static homogeneous 
and isotropic spatially flat cosmological solutions. This result has been ob¬ 
tained in Refs. |321I33]- 


5.1.2 First branch of bigravity 

From Eq. fl5.17l) we get a relation between the two Hubble constants and 
taking into account the Friedmann equations 05.101) . 05.lip we can derive the 
equation relating the two matter energy densities 


Hf = r-^H„ 


2 2m^ 

\ Blir) D 

Pg - pr pf + 

P Bo{r) 

K 

r 


where 



The hrst of Eqs. 05.20p gives 


pr^Pf + 3pf) + ^ (/x(Hi - 3rB2) + 3r^Bi) 
r [r{pg + 3p,) + ^ (r(Ro - 3Zlo) + 3/iZli)] ‘ 


(5.31) 

(5.32) 

(5.33) 


(5.34) 
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(5.35) 


Eq. fl4.24p aquires the following form 

r = (1 — ru)NHg, 
and a fixed point for r appears at 


_ 1 
u 


(5.36) 


It corresponds to proportional space-time metrics = r'^g^.v 


In general case for the expression 1 — wr we have the following formula 
(p, + 3p,) - pr2(p; + 3p/) + (ifi - Bo) 


ur = 


(P, + 3P,) + “ (Bo - 3B„ + iPSx) 


(5.37) 


In order to solve the system of Hamiltonian equations one can take as the 
initial data Pg(to) and Hgito). Then the constraint S 

An 


tj2 _ ^ I 

g QP3+ 


(5.38) 


allows to hnd Ag{to), and next, 'r(to) = ’"o can be found as one of the roots 
for the cubic equation 

Bo(ro) = 


Then constraint hli: 


Kr, 


Hj = r-^Hg, 


gives us Hf. The constraint TZ': 

- a'"’ +T- 


A/ = g,m 


2 S 

K, 


(5.39) 

(5.40) 

(5.41) 


provides us with pf. With account for Eq. (l5.37p and equations of state 
Pg — Pgip)y Pf — Pfip) it is possible to integrate the dynamical equations 

f = (1 — ru)NHg, 

Pg = -SNuHg{Pg+Pg), 

Pf = -^NHf{pf+pf), 

Nk, 


Hg = 


Hf = 


4 

NKf 


u{pg + Pg) - {1 - ur) - Di (r) 

hi 


Pf + Pf+i^- ur) 


2m? Di{r) 


K 
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where N{t) is an arbitrary monotonic fnnction responsible for a freedom of 
the time reparametrization. 

Besides the general case, when 1 — nr is given by Eq. fl5.37p . it is tempting 
to consider the case when both types of matter have the same eqnation of 
state 

Pg = wpg, Pf = wpf, (5.42) 

then Eq. fl5.37p is simplified, and is looking as follows 


1 — ur 


3(| + w) (/i^ - Bp) 
Bo-3Do+^-i^ + ^{l + 3w)pg' 


(5.43) 


5.1.3 Second branch of bigravity 

In the second case fl5.18l) the value of variable r is a constant: 

Zli(r) =/?i + 2/^2’^ += 0, -H- r = ^ . (5.44) 

P3 

Then from the kinematical equation fl4.24p it follows 


u = 


Hi 


(5.45) 


The last equation is not a constraint, in contrast to fl5.3ip . therefore the 
Friedmann equations (15.101) . (15.lip , and correspondingly, the dynamics for 
metrics decouple in this case. Dynamics for the two metrics are the 

same as in GR with the cosmological constants 


Ag = m^^Bo{r), 


Af = pm 


^Kg Bi{r) 


K m 


Eqs. (l5.23p . (I5.24p are as follows 


(5.46) 


TT Hhif 

Hf = . (Pf+Pf), 

(5.47) 


Hg — ^ [Pg + Pg) ; 

(5.48) 


and are the same as in GR (I4.42p . If we put ourselves into g-world, then the 
only artifact of the bigravity will be the hxed value of Ag. 
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The full set of equations consists of the algebraic ones 


r = 


= 

m = 


u = 


-/j2±^//3|-/3l/33 

/Sa 

^9 I -^9 

+ T’ 

+ T’ 


and the dynamical ones 


Y = 


Hf = 


NUKg 

~ 

Nkj 


{P9 + P 9 ) > 


(P/ + Pf) . 




Ul 

U) 


P9 


Pf 


3(P5+Ps) 2>{pf + Pf) 


= NH 


/• 


(5.49) 

(5.50) 

(5.51) 

(5.52) 

(5.53) 

(5.54) 

(5.55) 


5.2 g-matter without f-matter 

To obtain non-trivial results for the cosmology of bigravity which do not 
contradict the experimental evidence it is enough to study the case when 
f-matter is absent hrst considered in articles [2H1 ESI EO] . general, a solu¬ 
tion allows for an evolution starting from a matter dominated Universe and 
evolving to de Sitter geometry at late times. As there are many parameters 
/Si, i = 0,..., 4, and k/, the possible scenaria are diverse, including the cyclic 
Universe without the cosmological singularity. Even the elaborated statis¬ 
tical study presented, for example, in article [HI] does not alow to £x the 
values of free parameters. 

Let us suppose pf = 0 = p/. This case is discussed most often. As matter 
helds do not appear in Eq. fl5.15p for the secondary constraint there are 
the same two branches as in the previous subsection 

= ujHf — ^Hg = 0, (5.56) 

^2 = /3Y + 2/52^^^ +/SsY = 0. (5.57) 

The second branch gives nothing new, therefore we consider the hrst branch 
only. Suppose for simplicity k = Kg, then the constraints and dynamical 
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equations are as follows 


Pa ~ 

2m^ ( Hi r 

S V ^ / 

(5.58) 

Hi = 

Hi(r) 

P Q 

3 r 

(5.59) 

H, = 

N 

— [uK.gipg + Pg) + (1 - Ut) Hl(r)] , 

(5.60) 

Hf = 

4 (1 ^3 ’ 

(5.61) 

ur = 

(p» + 3p») + ^ ('F - B,) 

(5.62) 

(P, + 3p,) + (Bo - 3£>„ + Ml.) ■ 


Here we cannot freely specify both Pg, and Hg as the initial data. Only one 
of these quantities may be hxed, and then we are to solve quartic or cubic 
equation for r. Any variable, besides N, can be expressed as a function of r, 
and dynamics of r is given as follows 


r = (1 — ru)NHg, 

whereas N{t) is an arbitrary monotonic function. 

For an equation of state Pg = wpg we obtain 

, (1 + tt.) (gf - gp) 

1 _ i _ 

{l + w)i^-{l + w)Bo+{r + !^) Di’ 

then the dynamical equations aquire the following form 

(l + ») (A - B,) ^ 

T = I\ - - - 

Cs + w)l^-{l + w)Bo+{r + ‘^) Di’ 

H = (1 + ^^) + (^ + 3)rF>2] 

" 4 [l + yj)l^-{l + w)Bo+{^ + l)rD,' 


(5.63) 


(5.64) 


(5.65) 

(5.66) 


5.3 One matter minimally interacting with two met¬ 
rics 

Let us consider the scalar held having a minimal interaction with both f^i, 
and Pgi, metrics [SI ESI IMl ET] 

C4. = - «(■«) + \/=9 - «(■«) . ( 5 . 67 ) 
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In the homogeneous cosmology we suppose 0 = (j)(t) and 

TT^ 

^matter = p = / ^ + (m + r^)W(0) , (5.68) 

2^6 (- + r3) 

Then the Hamiltonian is as follows 


H 


H 


potential 


+ H 


kinetic 


H 


matter 


NTZ, 


(5.69) 


and the primary constraint-like equations are 


7^ 

5 




r^H] 


+ 




Kn 


2t?7^ 

+ p{u) + - {B,{r) + uBo{r)) 

K 


du ^ 


QH^ 2rn^ 

- ^-+p\u) + ^Bo{r) 


Kn 


K 


= 0 , 

(5.70) 


where 

P'(«) = (5-71) 

but both 77 and S nonlinearly depend on the auxiliary variable u. Then we 
should solve equation 5 = 0 for u and substitute the result into 77. Therefore 
this model contains only one (first class) constraint 77 = 0 responsible for the 
time reparametrization invariance, and so it is not free of the Boulware-Deser 
ghost, as also the model to be considered in Section [6l Here it is possible 
to specify more initial data than in the previous subsections, for example, p, 
Hg, and Hj. 


5.4 One matter minimally interacting with the effec¬ 
tive metric 

The new variant of interaction between the matter and bigravity is by means 
of the effective metric. When we do not restrict ourselves to mini-superspace 
the effective metric coupling leads to reappearance of the ghost as shown 
in Refs. [371 EHl EH EH SO]. Nevertheless it is claimed [371 EB] that there 
is a valid interesting region of applicability for this coupling at scales below 
the cut-off. Also, in a recent article [H] it is stated that the ghost can be 
excluded in a version of the terad formalism that is non-equivalent to the 
metric one. 
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The first analysis of cosmology in this model was given in article |12] . 
The effective metric is constrncted according to the formula 

+ /SVm., (5.72) 

in tetrad language this corresponds to exploiting the linear combination of 
the two tetrads: 

aE^ + (5.73) 

such a possibility for the special case a = f3 has been mentioned hrst in 
Ref. [25]. In minisuperspace the components of the effective metric will be 
denoted as follows 



Goo = Goi = 0, Gij = a^5ij = 'Gij, 

(5.74) 

where 

M = N{au + /3), a = + fiu. 

(5.75) 

There are 

the following standard relations: 



\/^ = N'y/Gp, \/^ = 

(5.76) 

where 

G = det{Gf,u), Ip = deti'ipij). 

(5.77) 

The primary constraints are 



S = Hg^ + + ^"''e"Ro(r), 

Kg K 

(5.78) 


Kf K 

(5.79) 


where it is supposed that the scalar held canonical variables (0, depend 
on time variable only and therefore 







(5.80) 


We introduce the following notation for the matter energy density 


P 





(5.81) 
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The minisuperspace Hamiltonian is as follows 

H = NiJZ' + uS), 


and the Poisson brackets are 


{«. H,] = 




{u,Hf] = 




12^2’ ^2a;2’ 

It is suitable to write Eqs. fl5.78p . fl5.79p in a different form 


(5.82) 


{4>,7i^} = l. (5.83) 


5 = e' 


7^' = oj^ 


6 


2m^ 


—Hg +cr(cr + /?r) p-\ - Boi^r) 


K. 


K 


6 o /3 ^ . n Bi{r) 


—Hf + —(a + /SrYp + 


Kf " K 


(5.84) 

(5.85) 


Demanding the fulhlment of equation iS = 0 in the course of evolution we 
obtain the secondary constraint, 


5 = {5,H} = Nn, 
where D has the following form 
Qm? 


-)■ 


D = 0, 


(5.86) 


D = {5,7^'} = ^ {uHf - iHg) [(3,e + , 


and we have introduced the notation for pressure 


vr; 


After calculating 




{D,5} = A ^ 0, 


(5.87) 


(5.88) 


(5.89) 


we see that D and S are 2nd class constraints. It is possible to introduce the 
Dirac brackets 


{F,G}n = {F,G}- 


{F,D}{5,G}-{F,5}{D,G} 

A 


(5.90) 


and consider TZ' as a hrst class contraint. The reduced Hamiltonian and 
equations of motion will be as follows 


H 


reduced 


— NTZ', F — {F, Hreduced} 


D- 


(5.91) 
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An equivalent way to derive the Hamiltonian equations is to exploit the 
Poisson brackets, but after their calculation insert a solution for the variable 
u determined from the equation 



n = {n,R} = N ({H, 7^'} + m{H, 5}) = 0, 

(5.92) 

i.e. 

u = -^{n,n'}. 

(5.93) 

As H 

is factorized 115.87lh compare Eq. (15.Kip') 



^ 2 5 

(5.94) 

there 

are two diherent solutions for it: 



Hi = ujHf — ^Hg = 0, 

(5.95) 

or 

^2 — + 2/320;^ + PsU P — 0. 

(5.96) 


Below we will consider both the two cases. 

With the Hamiltonian (15.82^ and the Poisson brackets (I5.83p we can de¬ 
rive the kinematical Hamiltonian equations, which, of course, are equivalent 
to the dehnitions of the Hubble constants (see fl4.23p to compare) and the 
scalar held momentum 

( O/T^ / O \ 

aff7 j “ 

Combining the kinematical and dynamical Hamiltonian equations for the 
scalar held (i.e. for the matter) 

0 = N(au + ns = —N(au -|- l3)a^U'{(j)) 

we obtain the energy conservation law 

^ = 0 , 

O'* 
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which can be written in the standard form 

p + 3-(p + p)=0. (5.100) 

a 

Dynamical Hamiltonian equations for the gravitational variables are the fol¬ 
lowing 


Hf = 


Nk 


f 


2m? 


f3{a + (3rf [p + p 


(3 + au 


+ 


H, = 


K 

NuKn 


(1 — ur) Di{r) 


a{a + (3r)^ ( p + 


p (3 + au 
w /d + 7 


2m^ .Di{r) 

-(1 — ur) 


K 


u 


(5.101) 

(5.102) 

(5.103) 

(5.104) 


(see fl5.23p . fl5.24p for comparison). Normally the matter equation of state is 
known also. 


5.4.1 First branch 

Here we consider the following solution of the constraint D = 0: 


Di = ojHf — ^Hg = 0, 

Then for the Hubble parameters 


Hf = 


Hn = 

" Nue 


Nu' 

the following relations are valid 




H = 


Hf = r-^Hg. 


N{au -|- /9)a’ 


Hg = rHf, Hf = H(^^ + /3y Hg = H{a + Pr). 


(5.105) 


(5.106) 


Constraints fl5.84p . fl5.85l) are equivalent to equations of the Friedmann form 




= 


Kjn 


a{a + (3r)^p 


2m^ 


Bn(r 


K 


(5.107) 

(5.108) 
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with a lot of new parameters a, /?, /Sq, ..., and a new variable r{t) involved. 
It is suitable to rewrite these equations as follows 


where 


Hi = 

+ T' 

(5.109) 

Hi = 


(5.110) 

Pf = 

+ P^ 

( 5 . 111 ) 

Pa = 

a{a - 1 - (3ryp, 

(5.112) 

= 

^Kf Bi{r) 
m — — 5 —, 

K 

(5.113) 

^9 ~ 

rn^—Boir). 

(5.114) 


From Eq. fl5.93p we obtain 


a - 3Z)2 + 3f,-irZ)i + + Hrf (pr-‘ + 

U = -y ^ , 

3F>i + - 3Do) + ^a(a + /3r)3 

(5.115) 

where 


Di — j3i + 2/3j+ir + /3i-^-2r‘^, 


Bi = Di + rDi 


j+i- 


(5.116) 


With account for Eqs. fl5.106p we can reformulate our problem to a couple of 
Friedmann-like equations for H 


= 


—a{a + I3r)p + —--—— 

6 3 K [a + pry 

tt2 n p n \ Bl 

H — —y {a +/3r) p + — --——rvw, 

6 r 3r k [a + pr)^ 

This system can be solved for the matter density 

,2 bMI _ BJr) 


2 m 

p = p r 




K (a -|- j3ry {p ^ra — (3) 


(5.117) 

(5.118) 


(5.119) 
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Vice versa it is possible to exclude p, and obtain 


^2 _ aBiir) - fiBoir) Kg 
3 (a + /3r)2 — (3) n' 

or we can express the cosmological term as a function of r: 


A 


= m 


Kg Bq 
K {a + (3rY 


(5.120) 


(5.121) 


When the unobservable variable r is given, we can calculte p, if, and A 
from Eqs. fl5.119p - fl5.12ip . Vice versa, if we take an initial value for p or if, 
after solving a quartic or cubic equation we can End a few solutions for r. 
The dynamical equation for r 


r = a(l 


ur) 


a + j3r 
(3 + au 


(5.122) 


is sufficient to predict its evolution, and therefore the evolution (forward or 
backward in time) of all observables. For the important factor 1 — nr we get 
the following expression 


1 — ur 


a{a + (3rY{p + 3p) (l - 
a{a + ftrf {p - 


M) ^ ^ - Bo) 

ar J _ K \ T _ ^ _ 

+ ^ (Eo - 3iio + ’ 


(5.123) 


Eg. (15.1 23 p may be simplihed if the matter has the equation of state p = wp\ 

3(l + n;)(^fi-i?o) 


1 — ur = 


- 3C„ + + (1^ - B„) (.^ + 


(5.124) 


r{t) may be non-monotonic and may change its behavior at those points 
where r = 1/u. 


5.4.2 Second branch 

The second branch 112 = 0 gives the following equation 

I3i(^ + 2l32ioj + l3oOj‘^ - = 0 ’ 


34 



















or, equivalently, 


2m?‘ Di{r) 

^ K a/3(a + /9r)2 

This equation gives a possibility to determine the pressure of matter p, and 
also the density p, through the matter equation of state, as a function of 
variable r. Then from constraints iS = 0, 7?.' = 0 we can hnd Hg and Hj. 
Evolution of these variables are provided by the Hamiltonian equations 

^9 = a{a + /gr)^(p + p), (5.125) 

Hf = + Prf{p + p). (5.126) 

The combination of Eqs. 05.1251) . 05.1261) gives the following relation 


4 

Hf 


—1 ^ 3 

p -ur. 


(5.127) 


6 Non-dRGT bigravity 


Let the potential in bigravity will be of non-dRGT type, for example, as 
taken from the Relativistic Theory of Gravitation (RTG) |15] : 


U 


- 1 




( 6 . 1 ) 


Then in 3 + 1 minisuperspace notations 

y/^ = Nue, ^f = Nu\ = - + ( 6 . 2 ) 


U 


and we obtain 




UJ 


3 1 — 


= NU, 


V 

w 


dU _ 
du 

U-u 


1 

2 

dU 

du 


2m2 




1 

2 





(6.3) 


(6.4) 

(6.5) 
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Therefore constraint-like eqnations are the following 




5 



2ni? 

H + 

- W = u' 


K 


2w? n 

n + 

V 




9 

m 

+ pf + — 
K 

2 

m I 





( 6 . 6 ) 

(^h.7) 


bnt in fact both Eqs. fl6.6p . fl6.7l) contain anxiliary variable u. Then we shonld 
solve one of these eqnations for u and snbstitnte the solntion into another 
eqnation. So, there is only one constraint in fact. It is a first class constraint 
arising as a resnlt of the time reparametrisation invariance having place in 
the bi-RTG. For this reason the variable N stays arbitrary. One can take, for 
example, pg, pf, Hg, and Hf as the initial data. We see that this theory has 
one more gravitational degree of freedom than the bigravity with the dRGT 
potential. As a resnlt, it is not snitable for the decription of the observable 
cosmology, and so never has been discnssed in the literatnre. 

The dynamical Hamiltonian equations here are as follows 


Hf = 

Nkj 

4 

H, = 

Nuk 

A 


9 


Pf+Pf + 


m? 1 — {urY 

K ur^ 


Pg+Pg- 


1 — {urY 

K, 


( 6 . 8 ) 

(6.9) 


6.1 Non-dRGT massive gravity 

When we take one metric, for example, as a background Minkowskian 
one, we put a; = 1, Hf = 0, p/ = 0, and the solution of Eq. fl6.6l) is 


1 


r'i 1^3 


( 6 . 10 ) 


After substitution of this into S we obtain the Friedmann equation in the 
RTG: 


tj2 _ '^gPg _ ^ ^g 
^ ~ 6 12^6 ^ 


(1 + - 35*). 


( 6 . 11 ) 


If we take n = Kg = 167rG, then it is possible to rewrite Eq. fl6.1ip as follows 


StiG 


H: = 

^ 3 


A 

-P, + 3 . 


( 6 . 12 ) 
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where 


(6.13) 


m 


nr ^ I 3 \ 

w)' 


We can take pg and Hg as the initial data. Then A(^) is determined from 
Eq. fl6.12p . and ^ conld be fonnd as a solntion of the bi-cnbic eqnation fl6.13p . 
In contrast to the massive gravity with dRGT potential this theory admits 
the homogeneons and isotropic dynamical Friedmann cosmology. It is due 
to the fact that the gravitational field has one additional degree of freedom 
here. But this degree of freedom is unfortunately a ghost. The effective 
cosmological constant is negative, and therefore any expansion should in 
future be transformed into contraction. The fate of the Universe is cyclic in 
this model. Of course, a quintessence matter held may provide an accelerated 
expansion for a hnite interval of time. 


7 Conclusion 

The cosmology of bigravity is an open problem that is under active study |6l 
[3 El EZ]. This article is written in order to demonstrate the beauty and 
power of the Hamiltonian formalism in this held of research. We would like 
to mention that it is a direct road to develop the quantum cosmology in 
bigravity. 

The author is grateful to V. A. Petrov, Yu. F. Pirogov, and Yu. M. 
Zinoviev for the interest to this work. 


8 Appendix 

The straightforward calculation of traces for matrix X starting from Fq. fl3.12p 
gives the following results 

TrX = -A + D, (8.1) 

TrX^ = - 2(HC') + TrT)^ (8.2) 

TrX^ = + 3A{BC) - (BDC) - TiD^. (8.3) 
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Then it is easy to estimate the function U and its derivatives 


V 

Vi 

w 


+ 


l3o& + l3ie ( X H—j- j + 1326 


+ Pze 


e + l 


+ 2 Trx^ + 


■^2 _ ^.^2 ^ 2 £!Lli££!£) 

e + l 

6{px^p) — Qx{px‘^p) — 3?/Tra:^ 
^ + 1 


+ 


-Ae(/p)* + 2(326 [{fx^pY - UpYx] + 

3/336 [ 2 x{fx^pY + Tra;^(/p)* - 2 (/(a:^)'^p)*] = -iiaCabPb, 

I3if + (3i6e + 2/3266 (x - ^ ^ + 

V £V + l)J 

6{px^p) — Qx{px‘^p) — 3zTrx^ 


+ 1^36 


— 3a;Tra;^ + 2Tra;^ + 


+ 1 


where 


Cab = 6 [Sab{(3i + 2(32X - 3/33Trx^) - 2xba{l32 + 3 / 330 ;) + 6(3sXbcXca] ■ (8.4) 


The following notations were used above 


^ab 

— f 

(8.5) 

X 

^aat 

(8.6) 

Trx^ 

^ab^bai 

(8.7) 

Trx^ 

^ab^bc^ca-) 

(8,8) 

Uab 

PaPc^cb') 

(8.9) 

y 

= {pxp) = PaXabPb = Vaa, 

(8.10) 

{px^p) 

Pa^ab^bcPa 

(8.11) 

(px^p) 

Pa^ab^bc^cdPd') 

(8.12) 

(ufp) 

= U^fiaVY 

(8.13) 

(ufx'^p) 

fia^baPb-) 

(8.14) 

ufix^p) 

fia^ba^cbPc' 

(8.15) 
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